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We search for the rare flavor-changing neutral-current decay B + — > K + vv in a data sample of 
82fb~ collected with the BABAR detector at the PEP-II B-factory. Signal events are selected 
by examining the properties of the system recoiling against either a reconstructed hadronic or 
semileptonic charged-B decay. Using these two independent samples we obtain a combined limit 
of B(B + — » K + vV) < 5.2 x 10 _J at the 90% confidence level. In addition, by selecting for pions 
rather than kaons, we obtain a limit of B(B + — » tt + v~v) < 1.0 x 10~ 4 using only the hadronic B 
reconstruction method. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Flavor-changing neutral-current transitions such as 
b —> svv and b — > dvv occur in the Standard Model (SM) 
via one-loop box or electroweak penguin diagrams with 
virtual heavy particles in the loops. Therefore they are 
expected to be highly suppressed. Because heavy non- 
SM particles could contribute additional loop diagrams, 
various new physics scenarios can potentially lead to sig- 
nificant enhancements in the observed rates 0,0. The- 
oretical uncertainties on b — > svv are much smaller than 
the corresponding b — > si + i~ modes due to the absence 
of a photonic penguin contribution and hadronic long dis- 
tance effects [3|. The SM B + — > K + vv branchingfrac- 
tion has been estimated to be (3.8±J;|) x 10" 6 jl 0, 
while the most stringent published experimental limit 
is B{B+ ->K+vv) < 2.4 x 1CT 4 at the 90% confi- 
dence level (C.L.) 0. There is additional suppres- 
sion of b — > dvv processes relative to b — > svv from 
the Cabibbo-Kobayashi-Maskawa matrix-element ratio 
\V td \ 2 /\V t3 \ 2 @. 

In this work we report the results of a search for 
the exclusive decay mode B + — > K + vv. By modify- 
ing the particle identification (PID) criteria used in the 
search, we additionally obtain a limit on the related de- 
cay B + — > tt + vv. Charge conjugate modes are included 
implicitly throughout this paper and all kinematic quan- 
tities are expressed in the CM frame (i.e. the T(45 l ) rest 
frame) unless otherwise specified. 

The data used in this analysis were collected with 
the BABAR detector at the PEP-II asymmetric-energy 
e + e~ storage ring. The results are based on a data sam- 
ple of 89 million BB events, corresponding to an inte- 
grated luminosity of 82fb _1 collected at the T(45) res- 
onance. An additional sample of 9.6 fb -1 was collected 
at a center-of-mass (CM) energy approximately 40 MeV 
below BB threshold. We used this sample to study con- 
tinuum events, e + e~ — > qq (q — u, d, s and c). Charged- 
particle tracking and dE/dx measurements for particle 
identification (PID) are provided by a five-layer double- 
sided silicon vertex tracker and a 40-layer drift chamber 
contained within the magnetic field of a 1.5 T super- 
conducting solenoid. A ring-imaging Cherenkov detec- 
tor (DIRC) provides charged K—ir separation of greater 
than 3<t over the momentum range of interest for this 
analysis, The energies of neutral particles are measured 
by an electromagnetic calorimeter (EMC) consisting of 
6580 CsI(Tl) crystals. The magnetic flux return of the 



solenoid is instrumented with resistive plate chambers in 
order to provide muon identification. A full BABAR de- 
tector Monte Carlo (MC) simulation based on GEANT4 [U 
is used to evaluate signal efficiencies and to identify and 
study background sources. 

The presence of two neutrinos in the final state pre- 
cludes the direct reconstruction of the B + — > K + vv sig- 
nal mode. Instead, the B~ meson from an Y(4S) — > 
B + B~ event is reconstructed in one of many semileptonic 
or hadronic decay modes, then all remaining charged and 
neutral particles in that event are examined under the as- 
sumption that they are attributable to the decay of the 
accompanying B. 

The B~ reconstruction proceeds by combining a D° 
candidate with either a single identified charged lepton 
or a combination A^~ ad of charged and neutral hadrons. 
The resulting semileptonic and hadronic charged-B sam- 
ples are referred to as B~j" and B^~ ad throughout this pa- 
per. The D° candidates are reconstructed by selecting 
combinations of identified pions and kaons that yield an 
invariant mass within approximately 3er of the expected 
D° mass in the modes K _ 7r + , K~ir + ir° and K~ir + ir + ir~ . 
For B^ ad reconstruction, D° — > A'g7r + 7r~ is also used. 

Photon candidates are obtained from EMC clusters 
with laboratory-frame energy greater than 30 MeV and 
no associated charged track. Photon pairs that combine 
to yield 77 invariant masses between 115 MeV/c 2 and 
150 MeV/c 2 and total energy greater than 200 MeV are 
considered to be ir° candidates. 

The B~, candidates are reconstructed by combining a 
D° candidate having a momentum p D o > 0.5 GeV/c with 
a lepton candidate of momentum pi > 1.35 GeV/c that 
satisfies either electron or muon identification criteria. 
The invariant mass, mui, of the D°£ candidate is re- 
quired to be greater than 3.0 GeV/c 2 . Under the as- 
sumption that the neutrino is the only missing particle, 
the cosine of the angle between the inferred direction of 
the reconstructed B and that of the lepton - D° combi- 
nation, described by the four vector (Eui^poe), is 



cos6>f 



2E h 



■ E 



Dl 



,D£ 



2\ PDi \- y/Ei 



beam 



(i) 



where mj is the nominal B meson mass and -Eboam 



and y/E\ 



beam 



are the expected B-meson energy 



and momentum, respectively. We use cos6b,d£ to dis- 
criminate against combinatorial backgrounds, for which 
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\cos6b.di\ can exceed unity. We retain events in the 
interval —2.5 <cos9b,d£< 1-1 in order to maintain effi- 
ciency for B~ — » D* g~D decays in which a 7r° or photon 
has not been reconstructed as part of the Di combina- 
tion. However events are vetoed if a charged tt consistent 
with B° — > D* + £~D is identified. If more than one D£ 
candidate is reconstructed in a given event, the candidate 
with the smallest \cos0B,De\ is retained. 

Reconstructed B^ ad decays are obtained by combin- reconstruction procedures result in raw yields of approx- 
ing a reconstructed D° candidate with a hadronic system imately 5800 B~^/fb~ and 2200 £^" ad /fb~ ■ Relative 
XT ad composed of up to five mesons (% ± ,K ± , and 7r°), in- systematic uncertainties of 4.5% (7%) are estimated for 



from B + B~ MC simulation against data. This pro- 
cedure compensates for differences in the B^ ad recon- 
struction efficiency in the low-multiplicity environment 
of B + — » K + vv events compared with the generic B + B~ 
environment. The Br, and ST reconstruction efficien- 
cies in MC are additionally validated by comparing the 
yield of events in which a B + — » D°l + v has been recon- 



structed in addition to the B had or B sl . 



TheB^andB^ad 



eluding up to two 7r° candidates. We define the kinematic 
variables m ES = yj ^bcam ~~ Ps and AE = E b - E beam 
where ps and Eb are the momentum and the energy 
of the B^ ad candidate. The system is selected by 

requiring that the resulting B^ ad candidate lies within 
—1.8 < AE < 0.6 GeV. If multiple B^ ad candidates are 
identified in an event, only the one with AE closest to 
zero is retained. The ibes distribution of reconstructed 
.Bhad candidates is shown in Fig. ^p. B^ ad candidates 
in the signal region, 5.272< m^s <5.288 GeV/c 2 , are 
used for the B + — > K + vv signal selection. Candidates 
in the sideband region, 5.225< mES <5.265 GeV/c 2 , are 
retained for background studies. 
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FIG. 1: (a) The D° mass distribution for D° -> K~n + decays 
used for B~[ reconstruction. Data are shown as points and the 
total background MC is shown as a solid histogram, (b) The 
mES distribution of B^ d events for data (points) and BB 
MC (solid histogram) . Continuum background has been sub- 
tracted from the on-resonance data using off-resonance data. 
The hatched histogram represents the estimated combinato- 
rial background from BB decays. 



Combinatorial backgrounds from continuum events are 
reduced in both the B~\ and B^ ad samples by requiring 
| cos6>t| < 0.8, where 8t is the angle between the thrust 
axes defined by the or B^ ad daughter particles, and 
by all other tracks and clusters in the event. Continuum 
events peak at | cos 9t\ = 1, while the distribution is ap- 
proximately flat for BB events. Backgrounds from QED 
processes are strongly suppressed by the B~ reconstruc- 
tion procedures and are negligible in this analysis. 

The B^ ad reconstruction efficiency for events con- 
taining a B + — > K + vV (signal) decay is determined 
from signal MC simulation after validating the yield 



the overall B~ x (-Bj~ ad ) yields. 

Events that contain a reconstructed B~ are examined 
for evidence of a B + — > K + vv decay. Tracks and EMC 
clusters not already utilized for the B~ reconstruction 
are assumed to be the daughters of the signal candidate 
B decay. Signal candidate events are required to possess 
exactly one additional track with charge opposite that of 
the reconstructed B~ . The track is additionally required 
to satisfy K PID criteria and to have momentum pk 
greater than 1.25 GeV/c. 

In addition to this track, B + — > K + vv events contain 
an average of approximately 200 MeV of EMC energy 
from hadronic shower fragments, photons from unrecon- 
structed D* — > D°j/tt q transitions in the B~ candi- 
date, and beam-related background photons. The total 
calorimeter energy attributed to the signal decay, i? cx tra, 
is computed by summing all EMC clusters that are not 
associated either with the decay daughters of the B~ or 
with the signal track. Signal events are required to have 
E cx tva < 250 MeV. The -E xtra distributions are shown in 
Fig. H for B s 7 



and B had events with one additional track 



that has been identified as a kaon. The B^ ad analysis 
additionally requires that there are six or fewer clusters 
contributing to -E C xtra, and that no pair of these clusters 
can be combined to form a tt° candidate. 

The total B + — > K + vv signal selection efficien- 
cies, including the B~ reconstruction, are esti- 
mated to be e K = (0.115 ± 0.009)% for B^ and 
e K = (0.055 ± 0.005)% for B had events. The quoted er- 
rors are the quadratic sum of statistical and systematic 
uncertainties. Theoretical uncertainties in the if ± en- 
ergy spectrum result in a 1.3% uncertainty on the signal 
efficiency. This uncertainty is evaluated by comparing 
the pk spectrum of B + — > K + vv MC events generated 
with a phase-space model with the models given in 0, 0] . 
Additional systematic uncertainties associated with the 
B + — > K+vv signal candidate efficiencies include the sin- 
gle track efficiency (1.3%), PID (2%) and EMC energy- 
modeling (3.8% for B A and 2.3% for B had ). The EMC 
energy modeling systematic is determined by evaluating 
the effect of varying the MC E extI3 , distribution within a 
range representing the observed level of agreement with 
data in events with a reconstructed B + — * D°l + v (for 
the B~, sample) and in samples containing two or three 
additional tracks (for the -B^~ ad sample). 
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Background events can arise either from B°B° or con- 
tinuum events in which the B~ candidate is constructed 
from a random combination of particles, or peaking back- 
ground events in which the accompanying B~ (or in the 
case of ST, at least the D°) has been correctly recon- 
structed. 

In the B~y analysis, purely combinatorial backgrounds 
are estimated by examining sideband regions of the re- 



constructed D° invariant mass distribution, 



de- 



fined by 3er < 



< 10<7 as is illustrated in 



Fig. for the D° — > K~n + mode. The sideband yields 
are scaled to the signal region under the assumption that 
the combinatorial component is flat throughout the D° 
mass distribution. This assumption has been validated 
using samples of events in which two or three tracks not 
associated with the B~ reconstruction are present. The 
total combinatorial background in the B~ x analysis is es- 
timated to be iV^- g > 3.4 ± 1.2. Although the peaking 
background prediction in the B~ x analysis have been stud- 
ied in MC and are shown in Figs and the peaking 
background in the final selection is not subtracted. 




by the MC statistical uncertainty (42%). Other system- 
atic errors include the overall B~ yield (7%), the remain- 
ing track multiplicity (5%), the particle mis-identification 
rates for the K± selection (6.3%), and the EMC -energy 
modeling (8%). The total peaking background in the 
Br, analysis is estimated to be 1.9 ± 0.9. The total 
(combinatorial+peaking) background in the B7 ad analy- 
sis is estimated to be = 3.9 ± 1.1 events. 
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FIG. 3: The pk distribution for (a) -Bj~ ad and (b) B^ events 
after applying the full B + — > K + uV selection except for the 
Pk > 1-25 GeV/c requirement. The data and background MC 
samples are represented by the points and solid histograms 
respectively. The dotted line indicates the expected signal 
distribution from MC 



The B + — ► K + vv branching fraction is calculated from 
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B{B+ -> K+uV) = 
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FIG. 2: The £ ex t ra distribution for B+ -> K' 
and B^ (right) events. Events are required to have a re- 
constructed B~ and exactly one additional track which has 
been identified as a kaon. No other signal selection cuts 
have been applied. The data and background MC samples 
are represented by the points with error bars and solid his- 
tograms, respectively. The dotted line indicates the expected 
B+^K+ vv signal distribution from MC 



In the Br, analysis, the combinatorial background 
can be reliably estimated by extrapolating the observed 
yields in the mss sideband region into the ttzes signal 
region, indicated in Fig. yielding 2.0 ± 0.7 events. 
The quoted uncertainty is dominated by the sideband 
data statistics, but includes also the uncertainty in the 
combinatorial background shape, which is estimated by 
varying the shape over a range of possible models. The 
peaking background in the BC ad analysis consists only of 
B + B~ events in which the -B^ad nas t> een correctly re- 
constructed, and is estimated directly from B + B~ MC 
simulation. MC yields are validated by direct comparison 
with data in samples of events in which the full signal se- 
lection is applied, except that either E extia > 0.5 GeV, or 
more than one track remains after the B~ reconstruction. 
Uncertainties in the peaking background are dominated 



where A^ bs is the total number of observed events in the 
signal region. N B ± = (88.9 ±1.0) x 10 6 is the estimated 
number of B mesons in the data sample and Ek is the 
total efficiency. A total of Af { bs = 6 (3) B+ -> K+vV 



candidate events are found in data in the B7 (Br,) 
analysis. The pk distributions for B + — > K + v~D signal 
events in the B~[ and Bj^ad analysis are shown in Fig. 
Branching fraction upper limits are computed using a 
modified frequentist approach, based on Ref. ||, which 
models systematic uncertainties using Gaussian distribu- 



tions. For both the B~ 



and B had searches. 



B+ -> K+i 



limits are set at the branching fraction value at which 
it is estimated that 90% of experiments would produce 
a yield that is greater than the number of signal events 
observed. Limits of B(B + — ► K + vv) s \ < 7.0 x 10~ 5 and 
B(B+ -> K + w) had < 6.7 x 10~ 5 are obtained for the 
B~y and Br &d searches respectively. Since the two tag 
B samples are statistically independent, we can com- 
bine the results of the two analyses to derive a limit of 
B(B+ -> K+vV) < 5.2 x 10~ 5 at the 90% C.L. 

We also report a limit on exclusive B + — > n + vv 
branching fraction using only the B had sample. The 
same methodology as for the B + — * K + vv search is ap- 
plied to the B + — > tt + vv search except that the single 
additional track is required not to satisfy either kaon 
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or electron PID criteria. The -E xtra and p„ distribu- 
tions for B + — > 7r + ^F are shown in Fig. 0] The over- 
all B + — > Tr + w selection efficiency is estimated to be 
c^r = (0.065 ± 0.006)%, where the quoted uncertainties 
include an estimated 2% PID uncertainty, and other 
contributions to the systematic uncertainty are similar 
to B + — » K + v~P. The peaking and non-peaking back- 
grounds are estimated to be 15.1±3.1 events and 9. Oil. 8 
events respectively, with similar systematic uncertain- 
ties to the B + — ► K + in7 analysis. The search selects 
N° hs = 21 candidates in data with an estimated total 
background of N^ 6 — 24.1 ± 3.6, resulting an upper limit 
of B(B+ -> 7r+w) had < 1.0 x 10" 4 at the 90% C.L.. 
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FIG. 4: The i5 cx tra (a) and p w (b) distributions for 
B + — » ■k + vV in the B^ ad sample. Events shown in the -Eextra 
distribution are required to have a reconstructed B~ and 
exactly one additional track satisfying the pion-selection re- 
quirements. The p n distribution has all signal-selection re- 
quirements applied other than the p w cut. The data and 
background MC samples are represented by the points and 
the solid histogram respectively. The dotted line indicates 
the expected signal distribution from MC. 

We see no evidence for a signal in either of the reported 
decay modes. The B(B + — > K + vV) limit reported here 
is approximately one order of magnitude above the SM 
prediction. It is the most stringent experimental limit 
reported to date. 
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